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Abstract 

We discuss the physical picture of thick vortices as the mechanism responsible for con- 
finement at arbitrarily weak coupling in SU{2) gauge theory. By introducing appro- 
priate variables on the lattice we distinguish between thin, thick and 'hybrid' vortices, 
the latter involving Z(2) monopole loop boundaries. We present numerical lattice sim- 
ulation results that demonstrate that the full SU(2) string tension at weak coupling 
arises from the presence of vortices linked to the Wilson loop. Conversely, excluding 
linked vortices eliminates the confining potential. The numerical results are stable un- 
der alternate choice of lattice action as well as a smoothing procedure which removes 
short distance fluctuations while preserving long distance physics. 
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1 Introduction 



Arguments for the presence of spread-out tubes of color-magnetic flux, thick 'vortices', 
being the essential feature responsible for maintaining confinement at arbitrarily weak 
coupling in SU (N) gauge theory were expounded some time ago [[J - . A key idea is 
that such extended structures cost very little action locally, and thus are not directly 
suppressed at large (3. By gradual variation of the gauge fields, they can disorder 
the vacuum over long scales. The infrared physical picture is, of course, independent 
of any ultraviolet cutoff details, but, as always with such nonperturbative questions, 
mathematically precise formulations have been possible only on the lattice. Thick 
vortices form closed extended structures which are topologically characterized, in the 
continuum extrapolation, by iri(SU(N)/Z(N)) = Z(N). 'Punctured' thick vortices, 
whose (small) 'hole' boundary is a Dirac monopole current loop, are also possible and 
survive at large (3 ||. (These Dirac monopoles are also classified by the non-trivial 
elements of 7ii(SU(N)/Z(N)).) In the SU(N) lattice gauge theory there also occur 
'thin' vortex excitations of the Z(N) part of the group. These are localized to one 
lattice spacing thickness, and hence are sensitive to the short distance details such 
as the precise choice of the plaquette action. Long thin vortices are very efficient at 
disordering the system at strong coupling, but are energetically heavily suppressed and 
become irrelevant at large f3. 

In this paper we first discuss in detail the various vortex excitations possible in 
the SU(N) lattice gauge theory (section 2). We treat explicitly the simplest N = 2 
case since no additional physical features appear in the general N extension which is 
straightforward. The proper distinction between thin and thick vortices and their inter- 
actions seems to have occasioned some confusion in the literature. A clean separation 
can be achieved by an exact rewriting of the SU(2) theory in terms of SU(2)/Z(2) ~ 
50(3) and Z{2) variables H, Jl|. We eschew any mathematical derivations that can 
be extracted from the literature, and give a detailed physical discussion of the various 
excitations and their interactions.^ In section 3 we examine the Wilson loop and its in- 
teraction with vortices. The expression for the Wilson loop operator in the SO(3)-Z(2) 
formalism shows that it is essentially a vortex counter. The fluctuation of the operator 
between positive and negative values is determined solely by the number (mod 2) of 
vortices linked with the loop. This then allows us to examine the vortex contribution 
to the string tension numerically (section 4). The string tension extracted from the 
full Wilson loops was compared to the same quantity extracted from the expectation 
of only the sign fluctuation counting the linking of the vortices. The computation was 
first performed with the Wilson action, and then also with a fixed point action in con- 
junction with a 'smoothing' procedure based on the renormalization group. The point 
of performing the comparison also under smoothing is that smoothing removes short 
distance fluctuations while preserving long distance physics. In particular, the string 
tension of the full Wilson loop remains unchanged under the smoothing procedure. A 

3 This is an extended version of the argument given in 0. 
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necessary test then of any claim concerning the long distance physics is that it remain 
invariant under the smoothing procedure. The numerical results demonstrate that the 
confining potential arises from the presence of the vortices linked to the loop: the full 
string tension is, remarkably, reproduced from the expectation of the vortex counting 
sign. Conversely, allowing no (mod 2) vortices to link with the loop eliminates the 
confining potential. Closely related results have been reported in [7j|. Our conclusions 
are presented in section 5. 

2 Vortices - thin, thick and hybrid 

For SU(N)/Z(N) gauge fields in the continuum, vortices are topologically classified by 
iii(SU(N)/Z(N)) = Z(N). This means that a vortex nontrivially linked to a Wilson 
loop (trace of the parallel transport matrix of the gauge field connection), taken in 
the fundamental representation of the covering group SU(N), contributes a factor 
z G Z(N), z ^ 1. A vortex forms a closed 2-dim surface in d = 4 (a loop in d = 3) 
so that it links with a Wilson loop C if it pierces once any surface bounded by C. 
Topologically, it is also possible to have gauge field configurations representing 'open' 
vortices (Dirac sheets). The boundary of an open 2-dim vortex sheet represents a 
monopole loop (monopole-antimonopole pair in 3-dim). These are Dirac monopoles, 
also classified by the non-trivial elements of tti(SU(N)/Z(N)) ||. 

Now in the continuum, where the gauge field is an element of the Lie group algebra, 
there is no local distinction between the pure SU (N)/Z(N) and SU (N) gauge theories. 
In the lattice formulation, in terms of group element bond variables, of course there 
is. The two differ by the dynamics of the additional Z(N) degrees of freedom present 
in the SU(N) case. Exciting these Z(N) degrees of freedom on a stack of plaquettes 
forming a 2-dim closed wall ( a closed loop in d = 3) gives a 'thin' Z(N) vortex. These 
are of course the vortices already present in a pure Z(N) lattice gauge theory (LGT). 
They are 'thin' because they necessarily have thickness of one lattice spacing. At small 
/3, they are very efficient at disordering the vacuum. At large j3, however, they are 
heavily suppressed by the SU (N) plaquette action, and get progressively frozen out as 
(3 increases. Correspondingly, the pure Z(N) LGT gets into a Higgs phase; whereas 
the distinction between the SU(N) and SU(N)/Z(N) LGT disappears, as it should, 
as the continuum limit is approached. Thus it is only the non-Abelian dynamics of the 
lattice analogs of the topological Z(N) = tti(SU(N)/Z(N)) vortices, which can, if at 
all, affect the large (3 long distance dynamics. 

The lattice literature contains many confused or incorrect statements due to failure 
to properly distinguish between the excitations of the Z(N) part versus the (lattice 
analogs of the) topological m(3U(N)/Z(N)) = Z(N) excitations of the SU(N)/Z(N) 
part of the SU(N) gauge group, and their respective energetics. A formalism that al- 
lows one to accomplish such a separation cleanly introduces separate SU(N)/Z(N) and 
Z(N) variables @,0j. It is important, of course, that this is done in a gauge-invariant 
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manner, and gives an exact rewriting of the partition function and all observables in 
terms of the new variables. 

From now on we restrict to SU(2), which is the actual case of our numerical simu- 
lations below. The extension to any iV is straightforward. Consider then the standard 
SU(2) theory partition function on a lattice A 

Z A = J n dU b exp (j2 P^ U )j > (!) 

where, as usual, we wrote U p = Ilbep Ub for the product of bond variables Ub around 
the plaquette p. 

We now introduce new Z{2) variables a p G {±1} residing on plaquettes. We write 
a c = Wp^cdp for the product of the o" p 's around the faces of the cube c. We also 
introduce the coset bond variables Ub G SU(2)/Z(2) ~ 50(3). The configuration 
space of the the SU(2) bond variables on the lattice A is split into equivalence classes, 
each class corresponding to one coset bond variable configuration {Ub} on A. Thus 
two ST/ (2) configurations {Ub} and {U' b } on A are representatives of the same coset 
configuration {Ub} if and only if one has U' b = Ub%, for some 75 G Z(2), for every 
bond b on A. Now given a coset configuration {Ub}, pick a representative {Ub} and let 
T] p = sign Up. Then the quantity 

Vc(U) = II Vp , (2) 
pec 

the product of r\ p around the faces of a cube c, depends, as indicated, only on the coset 
variables since it is invariant under Ub — > Ublb for lb G Z(2). In other words, it is 
independent of the representative used to compute it. 

Now one can show [I]], that ([!]) can be written in the form 

Z\ = I II dU b II d(J p II ^[Vco-c] exp(/3|trC/p|cr p ) . (3) 

J b V c 

In (|^), and what follows, the 'delta function' on Z(2) simply stands for 

5(r) = l[l+r], reZ(2) , (4) 

so that 8(t) = 1 for r = 1, 5(t) = for r = —1. Also, / da p (- ■ •) = =±1 (- • •) stands 
for 'integration' over the discrete Z(2) group. 

The crucial point is that the integrand in (|j) depends only on the f/'s since it 
is invariant under the local transformation Ub — * Ublb for arbitrary 7^ G Z(2). In 
particular, the action becomes the product of a Z(2) part and an S"0(3) part. The 
Z(2) part, which is given simply by the plaquette variable a p , determines the sign 
of the action. The 50(3) part is non-negative. In this connection note the relation 
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\tiU\ 2 = \xi/2(U)\ 2 = 1 + xi(U) = 1 + Xi(U), where Xi/2(U) and X i(U) denote the 
fundamental and adjoint representation characters of SU(2). Thus the [/-integration 
is in fact a [/-integration, i.e. in (^) one has 

J| dUb = Constant x J| dU b , (5) 

b b 

where dUb is the Haar measure over SO(3). 

Similarly, for the expectation of a Wilson loop one finds: 



W[C] = -Lj l\dU b txU[C] exp ^£ ^ tU p) 



(6) 



= -y- 1 n djj b n d °p n tiu i c ] ^ <?s ex P (^1^1^) . (7) 

J b P c 

In (0), U[C] = Yibec ^4 stands for the product of the U^s around the loop C, and we 
introduced the notations: 



Vs = II % , = II a p ' ( 8 ) 

pes pes 

where S is a surface bounded by the loop C, i.e. C = dS. It is easily seen that ([?]) does 
not depend on the choice of surface S. Furthermore, the quantity tr[/[C] rjs depends 
only on the SO (3) bond variables [/&, since it is invariant under £/& — > C47&, 7fe € ^(2). 
The Wilson loop operator is thus expressed as a product of an 50(3) and a Z(2) factor. 
Similarly, any other observable, such as a t'Hooft loop, or the electric and magnetic 
flux free-energies can be easily written in terms of the new variables. 

Equations and (|7|) then reexpress the SU(2) LGT as a coupled SO(3)-Z(2) 
theory. This rewriting is exact and gauge invariant. It is very convenient to evaluate all 
50(3) quantities in terms of SU(2) representatives, as in (|3|), (0).Q At the risk of being 
repetitive, let us point out again that, once the passage to the SO(3)-Z(2) formalism is 
made, the quantity r\ p = sign[/ p for a representative of a {Ub} configuration has nothing 
to do with the sign of the action in the 577(2) formalism ([[]), ©. This sign, as already 
noted, is supplanted by the Z{2) variable oy The representative-dependent r^'s can 
appear in physical quantities only in 50(3) representative-independent combinations, 
as e.g. (H). The resulting expressions fl3|), (0) of the SO(3)-Z(2) formulation have a 
physically rather transparent form exhibiting the presence and manner of coupling of 
the various possible topological excitations in the LGT ([!]). 

Consider a configuration where a v = 1 everywhere except on a stack of plaquettes 
forming a loop (Figure 0(a)) where a p = —1. This is a Z{2) vortex in 3 dimensions or a 
3-dimensional section of a vortex in 4 dimensions. In d — 4 there is an extra dimension 

4 They can always, of course, also be expressed directly in (a character expansion in) SO(3) (integer 
spin) representations, as indicated above for the action plaquette function. 
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Figure 1: Thin vortex: (a) d=3 or d=4 in [123]-section view; (b) d=4 in [130]-section 
3-dimensional view. Plaquettes protruding in x 2 -direction carry o p = —1. 

to move in, so by translation of the loop a vortex forms a 2-dimensional closed surface 
(Figure [j](b)). The short lines in Figure |I|(b) represent a set of bonds. The plaquette 
protruding in the ^-direction out of each of these bonds carries a p = — 1. A 3- 
dimensional [f-ivX] = [123] section gives then Figure |l|(a). These vortices, generated by 
the excitations of the Z(2) a p variables, are Z(2) 'thin' vortices, alluded to above . Note 
that, from (|3|), the action cost for exciting such a Z{2) vortex is directly proportional 
to the area of the vortex sheet. 

Consider next opening the thin vortex by breaking the loop of a = — 1 plaquettes 
in Figure 0(a) as depicted in Figure 0(a). The two cubes at the two ends necessarily 
satisfy a c — — 1. A cube with a c = —1 is the site of a Z{2) monopole, and Figure §(a) 
depicts a monopole - antimonopole pair joined by an open thin vortex, i.e. a string of 
a p = —1 plaquettes carrying the Z{2) flux. As just noted, there is a direct action cost 
associated with this cr-string. 

Now, because of the Z{2) 5-function constraint in the measure in (|3|), a cube with 
a c = — 1 must also have r\ c — — 1. A cube for which r] c = — 1 is the site of the 
lattice analog of a tt 1 (SO(3)) = Z(2) monopole. As pointed out above this statement 
depends only on the SO(3) coset configuration {Ub} on the lattice, i.e. the presence 
or absence of such a monopole on a given cube is a gauge-invariant feature of each 
SO (3) configuration. Any one representative {Ub} of an SO (3) {Ub} configuration 
with a monopole on a given cube will necessarily have a string of plaquettes, the Dirac 
string, beginning at the cube in question, on which rj p = —1. The string has to end 
at another monopole cube. Configurations with monopoles that contribute to the 
partition function then are of the form depicted in Figure |2](b) in 3 dimensions; heavy 
plaquettes carry o p = —1, light plaquettes carry rj p = —1. Note that the location and 
shape of the string (light plaquettes) depends on the choice of representative; it can 
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Figure 2: (d=3): (a) Z(2) monopole pair (cubes) joined by a-string (open thin 
vortex); (b) the complete configuration including the 77-string to form a 'hybrid' vortex 
(see text). 

differ for different representatives of the same {Ub} configuration since it may be moved 
around at will, as a Dirac string should, by letting Ub — > jbUb, % £ Z(2), i.e. change 
of representative. The ^-string is then "invisible" to the {Ub} configuration on A, and 
hence to the measure, and in particular to the action in (|3[), ([T]). There is no cost in 
action associated with the location of the Dirac 77-string. 




Figure 3: Hybrid vortex (d=4) in 3-dim [130]-section view. Plaquettes protruding in 
the remaining x 2 -direction out of light (heavy) bonds carry i] p = — 1 (cr p = —1); the 
boundary between the light and heavy sets is a Z(2) monopole current loop (see text). 

In d — 4 the monopoles, i.e. the set of cubes on which 1] C = —1, form closed 
monopole current loops reflecting magnetic current conservation. This follows directly 
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from the definition of the quantity T] c , eq. A string then sweeps out a Dirac sheet 

bounded by the corresponding monopole loop. So in d = 4 the configuration in Figure 
0(b) gives rise to a 'hybrid' vortex forming a closed 2-dimensional surface as shown in 
Figure |3|. Here again, a p = —1 on the plaquette protruding in the x 2 -direction out of 
every heavy bond, whereas r/j, = -1 on the plaquette protruding in the x 2 -direction 
out of every light bond shown. The plaquettes containing both a light and a heavy 
bond in their boundary are also shown in the figure. The cube protruding in the x 2 - 
direction out of each of these plaquettes then has a c = r\ c = —1. This set of cubes 
forms the monopole current loop (cp. footnote ||). A [123]-section view of Figure § 
gives then Figure 0(b). Such a hybrid vortex may be viewed as put together by joining 
an 'open' tti(SO(3)) vortex plaquette sheet and an open Z{2) vortex plaquette sheet 
along their respective monopole loop boundaries. These boundaries must coincide, as 
noted above, because open vortices as in Figure 0(a) cannot exist due to the constraint 
in the measure in (0). 

Closed Dirac plaquette sheets form the lattice analogs of tti(SO(3)) vortices (Figure 
£D, the set of plaquettes with r] p = — 1 being stacked over a closed 2-dimensional 
surface (a loop in d — 3). In the illustration of Figure | each plaquette carrying 
rjp = — 1 protrudes in the x 2 -direction out of the set of bonds shown distributed over a 2- 
dimensional surface placed in a 3-dimensional [013]-section. Again, the precise location 




Figure 4: Thick vortex closed Dirac sheet (d=4) in 3-dim [013]-section view; plaquettes 
protruding in the remaining x 2 -direction carry r\ p = — 1. 

and shape of this r\ v = — 1 set of plaquettes forming the Dirac sheet is irrelevant, it 
5 Indeed it follows from (|^) that rj c obeys the identity 

n ic=i , 

c e h 

where the product is over all cubes c forming the boundary of the elementary hypcrcube h. Geomet- 
rically, this means that the cubes on which rj c = — 1 form closed sets on the dual lattice. In d = 4, a 
cube is dual to a bond, so the cubes form closed loops of dual bonds. 
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being 'invisible' in the measure ([3]). What is relevant is only the coset configuration 
{Ub} describing the vortex; an SU(2) representative {Ub} of such a coset configuration 
will then contain somewhere on A a Dirac sheet, which may be moved around at 
will by a change of representatives. (Equivalently, the presence of the vortex can 
be characterized in terms of the Lf^s only - see below). Note that, since Z(2) (or 
generally Z(N)) flux is conserved only mod 2 (mod N), this implies that a vortex 
such as in Figure ^ is 'unstable' unless it is topologically nontrivial with respect to the 
lattice A or an externally introduced source in A, such as a Wilson loop. Indeed, by 
a change of representatives, the sheet may always be collapsed to a point annihilating 
the Z(2) flux, unless there is a topological obstruction. Thus on a lattice with periodic 
boundary conditions, i.e. the topology of the torus, a vortex as in Figure ^ may become 
topologically stable by wrapping completely around the lattice in the x 3 , x°-directions 
as shown schematically in Figure |5|(a). Here, short light lines represent a Dirac sheet of 
r] p = — 1 plaquettes, each in a [12]-plane, stacked along the x 3 , x°-directions around the 
periodic lattice. Such a topologically nontrivial closed sheet can be moved or distorted 
by a change representative, but not removed. Every representative of the relevant {Ub} 
vortex configuration has then an irremovable sheet of rj p — — 1 plaquettes signaling the 
trapped Z(2) flux of a (an odd number of) topologically nontrivial vortex (vortices). 
A characterization directly in terms of the Ub's is given by the quantity rjs defined as 
in (|D but with S now any closed topologically nontrivial surface winding around the 
lattice in the x , ^-directions (dashed line in Figure f|(a)). Then, clearly, rjs is only 
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Figure 5: (a) Topologically nontrivial vortex sheet (short lines) winding around periodic 
lattice (b) Vortex linked with Wilson loop C; 3-dimensional [130]-section view, with C 
lying in the [12]-plane. 



a function of the cosets {Ub}, and i]s(U) = — 1(+1) signifies the presence of an odd 
(even) number of vortices winding around the lattice in the x 3 , x°-directions normal 
to S. An analogous description applies to tti(SO(3)) vortices nontrivially linked with 
the Wilson loop (Figure 0(b)) discussed below. 
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As already noted, the Z(2) thin vortices are necessarily localized in thickness to 
one lattice spacing, and have a direct action cost proportional to their area (length 
in d = 3) along which a p = — 1. Long thin vortices are then heavily suppressed at 
large (3 as the a Z(2) variables are progressively frozen out. Indeed, it can be shown 
that the probability of exciting a p = — 1 on a plaquette is exponentially suppressed 
with large Only short thin vortices remain then with (exponentially in (3) small 
probability. This probability can actually depend on the choice of the lattice action.[| 
This dependence on the short distance structure is, of course, precisely a statement of 
the fact that thin vortices are thin. 

In contrast, the (lattice analog of the) tti(SO(3)) vortices are not necessarily lo- 
calized, and do not have a direct action cost proportional to their sheet area. Indeed, 
smooth {Ub} vortex configurations are easily constructed such that the local plaquette 
action cost can be made arbitrarily small by making the vortex sufficiently spread out 
|2|, IJ. Any SU(2) representative {Ub} of the {Ub} configuration for such a spread-out 
vortex will of course have a Dirac sheet plaquette set on which discussed 
above, while r] p = 1 everywhere else; but in such a manner that one still has \trU p \ ~ 1 
everywhere. The exact location of the Dirac sheet is in fact irrelevant since it can be 
moved by changing representative which does not affect |tr£7 p |. The point is, of course, 
that the action depends only on |tr£/ p |, i.e. {Ub}. These 'thick' vortices are thus not 
directly suppressed at large (3. In fact, long thick vortices winding around the lattice 
or a large Wilson loop can exist at arbitrarily weak coupling by being sufficiently thick 
in the directions transverse to their (topologically nontrivially linked) Dirac sheet. The 
same holds true for long hybrid vortices with a long thick vortex section, and a short 
(say, one plaquette long) thin vortex section appearing as a localized 'defect' incurring 
only a local cost in action |3|, 0]. These long hybrid vortices may be simply viewed 
as 'punctured' thick vortices, the size of the 'hole', where the a p Z(2) variables are 
excited, being suppressed, hence small at large (3. 

Such very long thick vortices (whole or punctured) can then have a very disordering 
long-distance effect at weak coupling as the bond variables Ub, over sufficiently large 
scales, vary smoothly over large parts of the SU(2)/Z(2) group with very little local ac- 
tion cost. Their presence appears in fact to be the necessary condition for confinement 
at weak coupling as we discuss in the next section. 

It is also interesting to view thick vortices in the ci-dimensional theory from a 
(d — 1) + 1-dimensional perspective by singling out the 'time' direction ||. The d-dim 
gauge theory may be viewed (in Kaluza-Klein (KK) fashion) as a (d — l)-dim gauge 
theory coupled to a Higgs field. A thick vortex may then be viewed within a (d— l)-dim 

6 This is proven nonperturbatively by 'chessboard estimates' [||. 

7 Thus the density of the thin Z(2) vortices, and/or Z{2) monopoles, may be enhanced or fur- 
ther suppressed by various short-distance modifications of the original lattice (Wilson) action in (|l]). 
Common modifications in the literature involve the addition of chemical potentials, the MP and the 
"positive plaquette" models, and models introducing new Z(2) degrees of freedom in addition to the 
aJs. 
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slice as a 'monopole'-'antimonopole' pair. Here the 'monopole' has two units (mod 2) 
of flux (and hence two 77-strings emanating from it), since it may be considered as put 
together out of two of our Z(2)-monopoles and is trivial under tti(SO(3)). Within 
a given (d — l)-dim slice, however, it may be characterized by using the 'Higgs' field 
of the KK dimensional reduction to define a homotopy n 2 (SO(3)/U(l)) group. These 
'monopoles' appear then as lattice analogs of the 't Hooft-Polyakov monopole in (d — 1) 
dimensions. In an appropriate gauge these correspond to the 'monopoles' of the Abelian 
projection. This is, of course, not a fundamental, gauge invariant description of the 
physical picture which is that of the iri(SO(3)) vortices. Still, it may be used as a basis 
for an approximate computational scheme for obtaining some estimate on the Wilson 
loop at weak coupling ||. Some numerical investigation of this picture has recently 
been reported in [0. 



3 The Wilson loop and Vortices 

Having identified the various types of vortices that occur in the SU(2) LGT, the ex- 
pression ([?]) for the Wilson loop is seen to have a rather transparent physical meaning. 
It makes explicit the interaction of the loop with vortices. Let us write (0) succinctly 

as 

WlC] = {trU[C]ns* S } somuz{2) , (9) 

where the expectation on the rhs is taken in the measure (^). As noted above, 
(^) decomposes the Wilson loop operator into a Z(2) part as and an 5*0(3) part 
tiU[C]r]s = \trU[C]\ sign (trC/[C] 775). It is crucial that the expectation (||) does not 
depend on the choice of the surface S spanning the loop. Then, for a given {{Ub}, {°~p}} 
configuration on the lattice: 

• If as = — 1 for every choice of the spanning surface S, a thin vortex, or an odd 
number of thin vortices, is nontrivially linked with the loop C (Figure [|(a)). 
Conversely, as = 1 for every S signifies an even number (including zero) of thin 
vortices linked with C . 

• If sign (trC/[C] 775) = —1 for every choice of the spanning surface S, a thick vortex, 
or an odd number of thick vortices, is nontrivially linked with the loop C (Figure 
|5](b)). Conversely, sign (trC/[C] 775) = 1 for every S signifies an even number 
(including zero) of thick vortices linked with C. 

• If neither as = —1, nor sign (tr[/[C] rjs) = —1 for every choice of S, but as sign (tr U[C] rjs) = 
— 1 for every choice of S, a hybrid vortex, or an odd number of hybrid vortices, 

is nontrivially linked with C . (Figure H(c)). Conversely, as sign (trC/[C] rjs) = 1 
for every S signifies that an even number (including zero) of hybrid vortices is 
linked with C. 
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Figure 6: Wilson loop linked with: (a) thin vortex, (b) thick vortex, (c) hybrid vortex; 
3-dimensional [123]-section view (cp. figure @(b)). 



Thus the passage to the variables Ub, a p exhibits the Wilson operator explicitly as 
a vortex counter. The fluctuation of the operator between positive and negative values 
is entirely due to the presence of vortices linking with the loop. The expectation of 
this fluctuation is what essentially determines then the behavior of the Wilson loop. 

The physically inessential |tr?7[C]|, which contributes only a perimeter effect, can 
in fact be eliminated from expression (|7|) by switching from the Wilson loop to the 
electric- flux free energy order parameter F e \ [JlO] defined on a lattice with periodic 



boundary conditions in all directions (A = T ). Expressed in our variables (|^), the 
electric-flux free-energy is given simply by: 



exp(-F el ) = {r ]s a s ) so(3)uz{2) , (10) 

where S is any 2-dimensional closed topologically nontrivial surface winding around 
the lattice in two given directions as in Figure 5. (The expectation does not depend on 
the specific choice of S.) By the above enumeration, (^) shows that the electric-flux 
free energy operator is nothing but a vortex counter for the various types of vortices 
winding around the lattice (in the directions perpendicular to S). 

Let us return to the consideration of the Wilson loop expectation (|9|). At large 
(3, long thin vortices, having large vortex sheet area A, disappear as they incur an 
action cost proportional to (3 A. Only a dilute gas of short thin vortices remains. For 
a large Wilson loop, short vortices can link with it only along the loop perimeter (as 
depicted for the thin vortex in Figure |5|(a)). At large /3, therefore, thin vortices can 
only contribute at most to a length-law piece in the expectation (0). This in fact can 
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be proven rigorously; it is equivalent to the statement that a pure Z(2) theory is in a 
Higgs phase at large (3. 

Long thick vortices, on the other hand, are not directly suppressed by the action 
at large j3. They may therefore link with a large loop anywhere over the area enclosed 
by the loop (as in Figure [1(b)). This may then lead to area-law behavior for the 
expectation, provided that the class of thick vortex configurations contributes at large 
P with a finite measure in the path integral sum. The same holds for long hybrid 
vortices, i.e. long thick vortices 'punctured' by a small monopole loop forming the 
boundary of a short thin vortex segment (Fig. ^](c)). Indeed, long Z{2) monopole 
loops are spanned by correspondingly large thin vortex sheets and suppressed at large 
f3\ but a dilute gas of short monopole loops survives at any finite (3. (The shortest 
possible loop is due to the excitation of a p = — 1 on a single plaquette p, forming a 
one-plaquette-long thin segment, and giving a Z(2) monopole loop consisting of the 
2{d — 2) cubes sharing this p on their boundary and hence having a c = —1.) In the 
absence of an artificial suppression of negative plaquettes (imposed, for example, by 
some modification of the action), this dilute gas of short Z(2) monopole loops can be 
used to tag hybrid vortices and estimate their contribution to the Wilson loop ||. 

One may modify the theory to exclude all Z(2) monopoles and hence all hybrid 
vortices by inserting in the measure (|3|) the constraint 

n %j • ( n ) 

c 

This is the MP model Confinement at large f3 in the MP model must then 

come from the thick vortices. Alternatively, one may instead eliminate all thick closed 
vortices linking with a given Wilson loop by inserting in the measure in the expectation 
(|) the constraint 

0[trtf[C]j7s] • (12) 

for any one particular surface S spanning the loop. Similarly, thick vortices winding 
around a periodic lattice may be excluded from the theory by inserting the constraint 

8[vs] , (13) 

for any particular closed topologically nontrivial surface S running through the lattice 
in two given directions (cp. Figure |5](a)): ([T3| ) eliminates all vortices winding in the 
(d-2) directions normal to S. In the presence of the constraints (|12"D or (|l3l), confining 

8 Note that the solution to the constraint (pi]), which is equivalent to requiring r) c = 1 on all cubes, 
is given by 

a P = II lh ' 

bep 

where jb are Z(2) bond variables, i.e. the Z(2) system in (|^) becomes exactly a (Wilson) Z(2) LGT. 
This is as one would expect: in the absence of Z(2) monopoles, only closed thin vortices are allowed 
as excitations of the <r p 's, which is the case in a pure Z(2) LGT. 
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behavior for (^), or (|TUD, respectively, at large (3 can then come only from the hybrid 
vortices. This is the approach taken in 0. 



Consider now inserting both the constraints (|TT1) and (|12[ ), resp. (|13|) in the mea- 
sure, thus eliminating both all hybrid vortices and all thick vortices winding through 
the Wilson loop, resp. the lattice. The form of the expectation (||), resp. fllOl) , now 
immediately suggests that confining behavior at large [3 is lost, i.e. that the presence 
of thick or hybrid vortices is the necessary condition for confinement to occur at weak 
coupling. In the case of the electric-flux free energy, eq. (|TUD, a mathematically rigor- 
ous proof of this fact was given in some time ago. The physical implications of this 
result for discussions of 'mechanisms of confinement' appears not to have been widely 
appreciated. It would clearly be important to have the corresponding proof for the 
case of the Wilson loop.[] Unfortunately, the proof in || does not immediately extend 
to the Wilson loop case. We will address this question elsewhere. 



4 Vortex contribution to heavy-quark potential - 
Numerical results 

As we saw in the previous section, the sign fluctuation of the Wilson loop operator 
is determined by its interaction with vortices. (In fact, in the case of the electric-flux 
free-energy, this interaction sign constitutes the entire operator.) This then allows one 
to directly examine the vortex contribution to the Wilson loop. We simply replace the 
value of the Wilson loop operator by its sign and consider the expectation 



E[C] = ( sign(trf/[C]) ) = ( sign ( trU[C] rj s ) a s ) so(3)uZ(2) . (14) 



The expectation ([[4]) is the vortex count expectation value as discussed above. In 
the following we wish to compare the string tension extracted from the full Wilson 
loop expectations (equations (0) and (0) ) with the string tension obtained from the 
expectation of the sign of the Wilson loops defined by equation ( |l"l"| ) . 

In all our measurements we extracted the heavy quark potential from timelike 
Wilson loops using the method and the code of Ref. p^| . We computed both on 
axis and off axis loops and the effective potential for different time extensions T was 
obtained as 

nSiT) = _ h TO (15 , 

v ' ; W(R,T) y 1 

9 Asymptotically large Wilson loops are physically essentially equivalent to the electric-flux free 
energy order parameter. In general, however, the string tension derived from the electric-flux free 
energy has only been rigorously shown to form a lower bound on the Wilson loop string tension [p"lf . 
Thus confining behavior for F e i implies confining behavior for the Wilson loop, but not, necessarily, 
the converse. In any case, as F e i is an order parameter that refers to the entire lattice, it is important 
to obtain the proof corresponding to the result in Q also for a large but finite Wilson loop as the 
lattice is taken to the thermodynamic limit. 
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In principle the heavy quark potential is the T — > oo limit of V(R, T). In the following 
we always display the effective potential for a time extent where it has already reached 
a good plateau. Typically with our values of the coupling this already happens at T 
equals a few lattice spacings. 

At first we used two ensembles of configurations generated with the Wilson action 
at (3 = 2.4 and 2.5 where the lattice spacing is a = 0.12 fm and 0.085 fm respectively. 
Our results are presented in Figs. (^) and (§). It is striking that the full Wilson loops 
and just their signs — the vortex expectations — give exactly the same heavy quark 
potential including the short-distance behavior and even the constant. We emphasize 
that we have not even shifted the two potentials by a constant, Figs. (|?p and (|8J) show 
the "raw" data without any further manipulation. 
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Figure 7: The heavy quark potential measured on an ensemble of 440 12 4 configurations 
generated at Wilson (3 = 2.4. Squares represent the potential obtained from Wilson 
loop averages, the octagons come from the sign averages. 
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Figure 8: The heavy quark potential measured on an ensemble of 129 16 4 configurations 
generated at Wilson (3 = 2.5. Squares represent the potential obtained from Wilson 
loop averages, the octagons come from the sign averages. 

The remarkable coincidence of the potentials computed in this way shows that 
the sign of the Wilson loop, i.e. the number of vortices (modulo 2) linking with it, 
contains all the important physics. The short- distance agreement of the potentials can 
be explained by noting that the sign expectation contains all the vortices including thick 
ones (77) and thin ones (a). The latter are important at short distances. Furthermore, 
as we saw in the previous section, thin vortices affect the perimeter-law term in the 
Wilson loop of any size. This in turn contributes to the constant term in the potential. 
The fact that even this constant is the same for the two potentials also shows that they 
contain the same contribution from thin vortices. 

At this point one could ask how robust this picture is, in particular how sensitive 
it is to the physically unimportant short distance details of the configurations. This 
can be checked either by modifying the action or by taking the Monte Carlo generated 
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configurations and performing some local smoothing on them which does not change 
the long-distance physical features. If the potential extracted from the sign of the 
Wilson loops is really equivalent to the full potential then their agreement at long 
distances should persist on the modified configurations. This is a very stringent test 
which has already been performed in the case of Abelian dominance. There it turned 
out that while on the original configurations the Abelian string tension agreed with 
the full SU(2) string tension to within 8%, after smoothing the difference increased 
to about 30% ||13||. Similar results have been obtained with cooling in Ref. WM. In 
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Figure 9: The heavy quark potential measured on an ensemble of 100 8 3 x 12 config- 
urations generated with a fixed point action (lattice spacing a = 0.14 fm.). Squares 
represent the potential obtained from Wilson loop averages, the octagons come from 
the sign averages. 

the present case at first we repeated the measurement of the full and "sign" potentials 
using the fixed point action of Ref. JTSJ at lattice spacing a = 0.14 fm. The results 



presented in Fig. (|9|) are very similar to the Wilson data; there is no measurable 
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difference between the potentials. We then performed one step of local smoothing on 
the same ensemble of configurations. This was done by the renormalization group 



based smoothing introduced in Ref. [fLq| . This local smoothing was designed to smooth 
only on the shortest distance scale, leaving all the long-distance physical features - 
most notably the string tension — unchanged. On the smoothed configurations the two 
potentials were measured again. Comparing the potentials obtained on the smoothed 
configurations (Fig. fllOD ) one can see that for distances R > 2 (in lattice units) they 
agree but for R < 2 the potential obtained from the signs is systematically below 
the full potential. This means that the smoothing destroyed a significant number of 
thin vortices and on this short distance scale thin vortices no longer dominate the 
the potential. On the other hand thicker vortices could not be destroyed by a local 
smoothing and the long-distance features are thus preserved. In this context we note 
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Figure 10: The heavy quark potential on the same ensemble as Fig. || but measured 
after one smoothing step. Squares represent the potential obtained from Wilson loop 
averages, the octagons come from the sign averages. 
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that exactly the same type of behavior would be expected from the positive plaquette 
model, in which the plaquettes are constrained to be non-negative. This constraint 
does not allow the formation of thin vortices but vortices thicker than one plaquette 
are not affected significantly. 
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Figure 11: The same as Fig. |9]but measured after 3 smoothing steps. Squares represent 
the potential obtained from Wilson loop averages, the octagons are obtained from 
Wilson loop sign averages. 

Finally we repeated the comparison of the full and the sign potential after an addi- 
tional two smoothing steps were performed (Fig. [II]). As a result of further smoothing 
the short distance disagreement of the potentials extended to a bit longer distances 
but the asymptotic string tension is not affected. This is consistent with our expec- 
tations that as more and more smoothing is performed, vortices of larger size are also 
destroyed. For a fixed number of smoothing steps, however, there is always a scale 
beyond which thick vortices remain intact. Beyond this scale one effectively has the 
same physical situation as before the smoothing. This may be viewed as being on a 
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fictitious coarser lattice with the lattice spacing set by this scale and with thin and 
thick vortices relative to this scale. Thus the vortex contribution to the asymptotic 
string tension is not affected by smoothing. 

5 Conclusions 

We presented a picture of the QCD vacuum by identifying the gauge field excitations 
that can disorder the system on large distance scales and can thus lead to confinement 
even at weak coupling. The relevant excitations are thick spread out center vortices 
that make the sign of large Wilson loops fluctuate considerably. The vortices are 
extended objects that cost very little in local action but have a long-range disordering 
effect. As opposed to thin vortices which gradually freeze out when the coupling is 
lowered, the thick vortices are expected to survive at arbitrarily weak couplings. 

We tested numerically how the vortices affect the Wilson loop expectations and the 
deduced heavy quark potential. In the SU(2) case vortices linking with the Wilson 
loop are responsible for the fluctuation of its sign. Therefore we compared the heavy- 
quark potential extracted from full Wilson loops with the potential extracted from the 
expectation of the sign of Wilson loops. The measurements were performed with the 
Wilson action at two different couplings as well as with a perfect action. In all three 
cases the two potentials completely agreed even for small distances. 

To check the universality of this picture we repeated the same test on an ensemble 
of locally smoothed configurations. The agreement of the long-distance part of the 
potentials persisted after the smoothing. This shows that all the relevant long-distance 
physical properties are encoded in the fluctuation of the sign of the Wilson loops which 
in turn is governed by the vortices linking with it. 
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